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ABSTRACT 

Topographic mapping of the earth, moon and planets can be accomplished with high 
resolution and accuracy using satellite laser altimeters. These systems employ 
nanosecond laser pulses and microradian beam divergences to achieve submeter vertical 
range resolution from orbital altitudes of several hundred kilometers. In this paper we 
develop detailed expressions for the range and pulse width measurement accuracies and 
use the results to evaluate the ranging performances of several satellite laser altimeters 
currently under development by NASA for launch during the next decade. Our analysis 
includes the effects of the target surface characteristics, spacecraft pointing jitter and 
waveform digitizer characteristics. The results show that ranging accuracy is critically 
dependent on the pointing accuracy and stability of the altimeter especially over high 
relief terrain where surface slopes are large. At typical orbital altitudes of several 
hundred kilometers, single-shot accuracies of a few centimeters can be achieved only 
when the pointing jitter is on the order of 10 (trad or less. 


i 



TABLE OF CONTENTS 


1. INTRODUCTION 

2. ALTIMETER RECEIVER MODEL 

3. GROUND TARGET MODEL AND LINK EQUATION 

4. MEAN PROPAGATION DELAY AND RMS PULSE WIDTH 

5. RANGE AND PULSE WIDTH ERRORS 


6. SURFACE SLOPE MEASUREMENTS 

7 BACKGROUND NOISE, THERMAL NOISE, QUANTIZATION NOISE 
AND SAMPLING EFFECTS 


8 EXPECTED PERFORMANCE CAPABILITIES OF SEVERAL 
SATELLITE LASER ALTIMETER SYSTEMS 


9. CONCLUSIONS 


REFERENCES 


Page 

..1 

..2 

...5 

...7 

.11 

.14 

.. 16 

..21 
..27 
... 28 


u 



1. INTRODUCTION 

Topographic mapping of the earth, moon and planets can be accomplished with high 
resolution and accuracy using satellite laser altimeters. These systems typically are pointed at 
nadir and function by measuring the round-trip propagation time of short laser pulses that are 
reflected from land, water or ice surfaces. The reflected laser pulses are distorted by surface 
variations within the footprint so that the shapes of the received waveforms also contain 
important information about surface roughness, slope and reflectivity. Meter or even centimeter 
vertical range resolution is possible from orbital altitudes of several hundred kilometers by 
employing gain switched solid state lasers with pulse lengths of 1-10 ns (0.3-3 m). Horizontal 
resolution is a function of the laser footprint diameter and pulse rate. Footprint diameters of 
several tens to several hundreds of meters and kilohertz pulse rates are readily achieved with 

current laser technology. 

Laser altimetry, especially from aircraft, has a long history of development and application. 
The first spacebome altimeters were ruby laser systems flown on the APOLLO 15, 16 and 17 
missions to the moon [Kaula et al., 1974]. During the past two decades numerous airborne 
system have been developed to measure geologic features such as volcanoes [Garvin et al., 

1991], ice sheet topography [Bufton et al., 1982] and sea state [Abshire and McGarry, 1987; Tsai 
and Gardner, 1982]. Bufton [1989] provides an excellent overview of the laser altimetry field 
while Bufton et al. [1991] describe recent technology advances that are yielding significant 
improvements in the accuracy and resolution of airborne systems. 

Several satellite altimeters are now under development at NASA, to study the earth, moon 
and Mars. These include the Geoscience Laser Ranging System (GLRS) which is being 
developed for the Earth Observing System [Cohen et al., 1987], the Mars Observer Laser 
Altimeter (MOLA) which will be launched in 1992 [Abshire et al., 1991], the Lidar In-Space 
Technology Experiment (LITE) planned for launch on the shuttle in late 1993 [Couch et al., 
1991], the Lunar Observer Laser Altimeter (LOLA) and the Topographic Mapping Laser 
Altimeter (TMLA) which is being developed as a small dedicated Earth Probe Class Satellite 
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[Harding et al„ 1991] . All of these instruments will to some degree, record information about 
the reflected pulse shape as well as the standard time-of-flight data. The use of laser waveform 
data in topographic applications has not been explored thoroughly. Gardner [1982] and Tsai and 
Gardner [1982] have studied the effects of the surface profile and sea state on the shapes and 
widths of ground and ocean reflected laser pulses. In this paper we extend their results to 
analyze the performance of satellite laser ranging systems. We develop detailed expressions for 
the range and pulse width measurement accuracies and use these results to evaluate the expected 
performances of the LITE, TMLA, GLRS, MOLA and LOLA altimeters. Our analysis includes 

the effects of the target surface characteristics, spacecraft pointing jitter, and waveform digitizer 
characteristics. 

2. ALTIMETER RECEIVER MODEL 

The geometry of the laser altimeter and ground target is illustrated in Fig. 1. We assume the 
ground target is diffuse and there are no specular reflections. Therefore, the analysis applies to 
reflections from land or snow but does not apply to reflections from water or ice. The detector 
model is illustrated in Fig. 2. A portion of the reflected optical pulse p(t) is collected by the 
receiving telescope and focused onto a detector. The detector is modeled as an ideal detector 
followed by a linear filter with an impulse response given by h(t). The detector output current 
i(t) is sampled and quantized by an A/D convener and then processed by a range computer which 
computes the signal level, pulse propagation delay and pulse width. 

The detected signal is contaminated by several types of noise which effect the accuracies of 
the computed pulse parameters. The optical signal is affected by both photon noise and speckle 
noise while the detector output will be contaminated by thermal noise and perhaps multiplication 
noise if the detector is an avalanche photodiode. Finally, the quantized signal at the output of the 
A/D converter will be contaminated by quantization noise. The effects of all these noise sources 
on the performance of the laser altimeter is considered in the following sections. 
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Receiving Transmitting 



Figure 1. Geometry of the laser altimeter and ground target for nadir pointing. 



(N,T p , O p ) 


Figure 2. Detector model for the laser altimeter. 

Ranging performance is influenced to a great extent by pulse shape, signal level and the 
algorithms used to compute the pulse delay. Unfortunately, minimum mean-square error (mmse) 
estimation techniques cannot be used with laser altimeters because the received pulse shape 
depends upon the surface profile and cannot be predicted a priori. The most effective technique 
for estimating range when the pulse shape is unknown is to compute the centroid of the received 
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pulse. The performance of this estimator is proportional to the rms pulse width and is inversely 
proportional to the square root of the signal photon count. In this paper we are primarily 
interested in deriving theoretical expressions for the pulse centroid time T p , the rms pulse width 
a p and the pulse photon count N. These three parameters arc functions of the temporal moments 
of the pulse and are defined as follows 

N=fp ( t )d t (1) 

OO 

T P = ^/tp(t)dt 2 

0 y ’ 

oo 

= n / (* ~Tp)2 p(t) dt (3 ) 

where p(t) is the expected received pulse shape, i.e. photon count rate (photons/sec). Of course, 
the receiver computes estimates of N, T p and c p 2 from measurements of the detector output 
current i(t). 

If we assume that the detector bandwidth is much larger than the signal bandwidth, then h(t) 
can be modeled as an ideal integrator 


e/At Itl < At/2 
0 otherwise 


(4) 


where e is the electron charge and At is the sampling interval or range bin length of the A/D 
converter. The expected output current of the detector is 

t+Al/2 

E[i(t)] = G p(t)*h(t) = G j p(t) dx ^ Ge p(t) 

t— At/2 

G = detector gain. 
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Since E(i) is proportional to p(t), the parameters N, T p and o p 2 can be computed from (1M3) 
with p(t) replaced by i(t) 

oo 

n=4t J i(t)dt 
06 0 

OO 

T P = GeN J tl(t ^ dt 

0 P 2 = ok J(t-Tp) 2 i(t)dt. 


In subsequent sections we derive theoretical expressions for the means and variances of these 
parameters. 


3. GROUND TARGET MODEL AND LINK EQUATION 

Since the ground target is a diffuse reflector, the statistics of the reflected signal are identical 
to statistics for fully developed speckle. The signal amplitude is a circular complex Gaussian 
process and the signal intensity obeys negative exponential statistics. The mean and variance of 
the total detected photon count are given by [Gardner, 1982] 


<N> = E 


Ge 


J i(t) dt 


hv nr 


(9) 


Var (N) = F<N> + <N> / K s 


GO) 


K s = TtA 


'2tan0T N ? 


(ID 
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where 


F = excess noise factor of the detector 
T a = 1-way atmospheric transmittance 
t| = detector quantum efficiency 

E = total transmitted pulse energy (Joules) 
hv = photon energy (Joules) 

A = receiver aperture area (m 2 ) 
r = z/cos0 = ground target range (m) 
z = altimeter altitude (m) 

<() = nadir angle of laser beam 
P(<j)) = diffuse reflectivity of the target 
0T = laser beam divergence, HW @ e-^ 2 (rad) 

X = laser pulse wavelength (m). 

K s is the ratio of the receiver area to the speckle correlation area and is usually called the speckle 
signal-to-noise ratio. The speckle correlation area is inversely proportional to the area of the 
laser footprint. Typically, K s is on the order of 10 4 - 10$ and <N> is on the order of 103 - lO 4 
The excess noise factor F is equal 1 for photomultiplier tube (PMT) detectors and between 3 
and 7 for avalanche photodetectors (APD). In general, the altimeter will not always be pointed at 
nadir. The system geometry for arbitrary nadir angles is illustrated in Fig. 3. The coordinate 
system is defined by the optical axis of the altimeter and the line drawn from the altimeter to the 
center of the earth. The nadir angle <> is the angle between the optical axis and center line. The 
center line is defined as the z axis. The 2-D surface profile is described by the function $( C ) 
where £ - (x,y) is the horizontal position vector on the ground which is normal to the z-axis. 

The surface profile within the laser footprint is modeled as 

^(c) = + £ • a + a^(ji) (l2 \ 
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where Z# is a constant offset, £ is the mean surface slope within the footprint and AS, is the 
surface roughness. The model given by (12) differentiates the linear variation in surface height, 
+ from quasi-random surface roughness variations AS,. By defining the origin of the 

coordinate system as the intersection point of the center line and the earth’s surface we have 
^ o = 0 and 


^(x,y) = x tanSii + y tanSj. + A^(x,y) 


(13) 


where Su is the surface slope parallel to the nadir direction and S± is the surface slope nomtal to 
the nadir direction. Similarly, we define Afll and Adi as the orthogonal components of the 
pointing jitter, parallel and normal to the nadir direction. 


4. MEAN PROPAGATION DELAY AND RMS PULSE WIDTH 

Expressions for the mean pulse propagation delay and rms width can be derived using the 
same assumptions and the approach employed by Gardner [1982]. If we assume the laser beam 
cross-section is Gaussian (TEM^ mode), then for the geometry illustrated in Fig. 3, the mean 

pulse delay is 

2z(l + tan 2 9j) [ 1 + tan 2 (<|) + A<tm ) + tan 2 A<|>x] 1/2 
E(T p I A<|>ii, A<t>i) - c cos( j ) [i _ tan(<j> + A<t>ii) tanSii - tanA<])j. tanSi] 


cos(<|) + Su) 
cosSh 


(14) 


and 


E(T p ) _ + tan £l + (1 +2 tan 2 (<|) + Su)) 

p c cos<]> 


V ar(A<J)||) 


+ COS 2 <() 


fl + 


V 



(15) 
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S(x,y) — x tan S|.+ y tan S|| + A^(x,y) 


Figure 3. Generalized geometry for the laser altimeter and ground target $(x,y) is the surface height, 
where 

c(<J>) = effective velocity of light (m/s) 

S|| = surface slope parallel to nadir direction 

Si = surface slope normal to nadir direction 
A0 h = pointing error parallel to nadir direction 
= pointing error normal to nadir direction. 
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The effective velocity of light depends on the laser wavelength, nadir angle and surface 
meteorological conditions. The additional roundtrip propagation delay introduced by the earth s 
atmosphere is approximately 16 ns (4.7 m) at nadir and is proportional to l/cos<|> [Gardner, 

1977]. 

The expected propagation delay given by (15) is composed of three terms. The 2z/(c cos<t» 
term represents the roundtrip propagation delay along the center of the laser beam. The term 2z 
tan 2 9r/(c cos<{>) is the additional delay resulting from the phase front curvature of the diverging 
laser beam. Finally, the terms involving Var(A<|>ii) and Var(A<t>i) are biases caused by pointing 
jitter. Equation (15) was derived by expanding Eq. (14) for the conditional mean delay in a 
power series in terms of A«Dj. and A<t>u, retaining the terms out to second order in both A0± and 
A(j>|| and by assuming the orthogonal components of the pointing jitter are statistically 

independent. 

For nadir pointing the range biases, <i.e„ cT,/2 ) introduced by beam curvature and pointing 
jitter are approximately 

Beam curvature range bias - z tan 2 0j ^ 

Pointing jitter range bias ~ z Var(A<])) . ^ ‘ 

The pointing jitter bias dominates whenever the rms jitter exceeds the beam divergence. For an 
altimeter altitude of 1000 km and an rms pointing jitter of 10 [trad, the range bias caused by 
pointing jitter is only 0.1 mm. For a beam divergence of 64 prad HW @ e~ 1/2 (150 prad 
FWHM), the range bias caused by beam curvature is 4 mm. 

An expression for the mean-square pulse width can also be derived for the geometry 
illustrated in Fig. 3 by following the approach of Gardner [1982] 
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Table 1 Typical Values for Pulse Broadening Contributions 


Source 

Simplified Formula 

Value* 

System effects 

§«Jt 2 + 0|,2)W 

~ 1 m@ Oh^ot = 5 ns 

Surface roughness 

Std(A^) 

0.01 — 10 m 

Beam curvature 

Z0j2 

<0.01 m@ z= 1000km 

Nadir angle effects 

z0j tan<j> 

1.1 m@ <J> = r 
5.5 m<5> <|> = 5* 
12 m @ <J> = 10“ 

Surface slope 

z0t tanS 

1.1 m@ S = 1* 
5.5 m @ S = 5“ 
12 m @ S = 10° 
37 m @ S = 30“ 


*AU values were computed assuming z = 1000 km and 0r = 64 prad (i.e., 150 prad FWHM) 


E(a p 2)= ( 0e 2 + 0h 2 )+ lVH(Mlco£Sj 

c 2 cos 2 (<t> + S||) 

(system effects) (surface roughness effects) 


4z 2 tan 2 0r 
c 2 cos 2 ^ 


tan 2 0x + tan 2 (<> + S||) + ffl 125 -!- cos 2 Sh 

COs2(<t> + S||) 


(beam curvature effects) (nadir angle and surface slope effects) 


(18) 


where 


Gi = rms laser pulse width (s) 

Oh = rms width of receiver impulse response (s). 


The system contributions to the width of the received pulse include the widths of the transmitted 
pulse and receiver impulse response. For wide bandwidth receivers which employ A/D 
conveners to sample the detected signal, o h = At/Vl2 where At is the sampling period. The 
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received pulse is also broadened by the surface roughness and slope of the ground target and by 

c 

the beam curvature. Typical values for the pulse broadening contributions (i.e., ^ °p) are listed 
in Table 1. Surface roughness, surface slope and nadir angle effects can be quite significant and 
may contribute several meters to several tens of meters to the pulse width. 


5. RANGE AND PULSE WIDTH ERRORS 

Photon noise, speckle noise and pointing jitter all contribute errors to the range and pulse 
width measurements. Again by generalizing the results of Gardner [1982], we obtain the 
following lengthy expression for the variance of the pulse propagation delay. 

F(q<t 2 + <7 h 2 ) / F l \ 4Var(A^) cos 2 S|| 

Var(Tp) <N> (^<N> K S J c 2 cos 2 (<)> + S||) 

(system effects) (surface roughness effects) 


/ F 1 \ 4z 2 tan 2 6>T 
1^<N > 2K s J c 2 cos 2 <J> 


tan 2 0j + tan 2 (<}> + S||) + 


tan 2 Sj_ cos 2 S||" 


4z 2 (1 + tan 2 Q-r) 2 
c 2 cos 2 <> 


cos 2 (<l) + Sh) 

(beam curvature effects) (nadir angle and surface slope effects) 

tan 2 Si cos 2 Sii cos 2 0 
cos 2 (<|> + S||) 

(pointing jitter effects) 


tan 2 (<t> + S||) Var(A<|>ii) +- 


Var(A<{)i) 


(19) 


The photon noise contributions to the range error variance are inversely proportional to <N>, the 
expected number of detected photons per received pulse. The speckle noise contributions are 
inversely proportional to K s , the ratio of the receiver aperture area to the speckle correlation area. 
Some of the speckle noise effects can be reduced by increasing K s either by increasing the 
aperture area or the laser divergence (see Eq. (11)). However, some of the speckle noise terms 
are proportional to tan 2 0x/Ks and tan^0r/Ks* These terms are either independent of 0 t or 
proportional to tan 2 0j. It is not difficult to show by taking the derivative of (19) with respect to 
tan0T that speckle errors are minimized when 
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Table 2 Range Measurement Errors | Std(T p ) 


Source 

Simplified Formula 

Value* 

Surface roughness 

Std(A$) 

V<N> 

negligible to ~ 0.6 m 

Nadir angle 

z0t tan<|> 

~ 0.2 m @ <[> = 5° 

V<N> 

Surface slope 

zQj tanS 
V<N> 

~ 0.04 m @ S = 1° 
~ 0.2 m (a) S = 5° 
~ 0.4 m @ S = 10° 
~ 1 m @ S = 30° 

Pointing jitter 

zStd(A<)>) tanS 

~ 0.18 m @ S = 1° & A<|> = 10 pxad 

~ 0.8 m @ S = 5" & A<|> = 10 (xrad 

~1.8 m @ S = 10° & A<(> = 10 (irad 

- 6 m @ S = 30' & A<|> = 10 prad 


tanQ T --l2 C ° SS " cos4> Std < A 5> 
cos(<)> + Sn) z 


( 20 ) 


The laser divergence is typically much larger than the value given by (20). 

The major ranging error effects | Std(Tp) are listed in Table 2 for the case where F = 1. 
Typical values for these errors were computed by assuming <N> ~ 10^ and K s ~ 10^. The 
surface roughness, nadir angle and surface slope contributions arise because the performance of 
the centroid estimator is proportional to the received pulse width. Pulse broadening is greatest 
over rough terrain or when the nadir angle or surface slopes are large. The pointing jitter 
contribution arises from the bias introduced into the propagation delay when there is a pointing 
error. Notice from the simplified formulas included in Table 2 that pointing jitter and surface 
slope contributions are comparable when 

Std(A<j>) = 9r / VN\ (2 i 
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£ 

Table 3 Pulse Width Measurement Errors ^Stdfap) 


Source 

Simplified Formula 

Value* 

Surface roughness 

Std(A^) 

V2<N> 

negligible to ~ 0.4 m 

Surface slope 

z0j tanS 
V2<N> 

~0.2m@S = 5‘ 
-0.3 m@ S = 10” 
- 0.8 m @ S = 30' 

Nadir angle 

z0t tan<(> 
V2<N> 

~0.2m@ <t> = 5* 


♦All values were computed assuming z = 1000 km, 0 T = 64 prad (i.e., 150 prad FWHM), <N> - 10 3 and 
F= 1. 


For <N> = 10 3 and 0x = 64 prad, we see from (21) that pointing jitter contributions dominate 
when Std(A<|)) > 2 prad. Pointing jitter is usually the dominant source of range error especially 
over terrain with large surface slopes. The extreme sensitivity to pointing errors is due to the 
small laser divergence. Terrain slopes averaged over the footprint can be quite large when the 
divergence is small. When divergence and hence the footprint are large, the average surface 
slope within the footprint is usually small, however the roughness can be large in this situation. 
Depending on the statistical characteristics of the surface profile, there is an optimum divergence 
angle which minimizes range error by balancing the surface roughness and the pointing 
jitter/surface slope contributions. 

If we neglect the effects of speckle, the variance of the rms pulse width measurement is given by 
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v^'-ESSl 2 ) F(g t 2 + °li 2 ) F 2Var(A£) cos 2 Si 
P 2<N> - 2<N> + <N> C 2 C0S 2 W + S||) 

(system effects) (surface roughness effects) 

F 2z2 tan 2 0T 


<N> 


C z CO 


~rr^ [ tan 2 0 T + tan(<> + Sh) + — - S± cos2s »l 
L cos2(<t» + S,|) J 

(nadir angle and surface slope effects) 


(beam curvature effects) 


( 22 ) 


The performance of the pulse width estimator is proportional to the received pulse width and so 
the major contributions to the error arise from surface slope and nadir angle effects. Typical 
values for the pulse width errors | Std(a p ) are listed in Table 3. 


6. SURFACE SLOPE MEASUREMENTS 


For some geophysical applications, the accurate determination of surface slope is of 
considerable interest. The total slope can be estimated from the measured pulse width. Consider 
the diagram in Fig. 4a for nadir pointing. The laser footprint diameter is approximately 
2z tan0 T tanS. The surface slope is approximately equal to the pulse width divided by the 
footprint diameter. From Eq. (18) with the nadir angle equal to zero, we have 


- ian2 S + tan^ „ VggB. , ° 2 «* 2 + °h 2 > 

4z 2 tan 2 0 T z 2 tan 2 0r 4z 2 tan 2 0 T 


(23) 


The last 3 terms in Eq. (23) are bias terms. The biases related to the laser divergence, laser pulse 
width and receiver impulse response can be calculated from known system parameters and 
subtracted from the right-hand side of Eq. (23). The bias related to surface roughness cannot be 
determined a priori but can be minimized by increasing the beam divergence. Of course, the 
range bias and error increase with increasing Qj (see Eqs. (19) and (20)). When z = 1000 km, 0 T 
= 64 prad and Std(A^) = 10 m, the slope bias is approximately 9*. The slope bias is less than 1° 

when the rms roughness is 1 m. The variance of the slope estimate can be calculated by using 
Eq. (22). The result for <f) = 0 is 
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Laser Beam 



Figure 4. Target and altimeter geometries for surface slope measurements: (a) single pulse technique and (b) 
double pulse technique. 


Var(tanS) = 


F 

2<N> 


tan 2 S + tan 2 6r + 


Var(A^) 
z 2 tan 2 0r. 


(24) 


If surface roughness effects arc negligible, the slope error reduces to 


Std(tanS) ~ tanS 


-/m 

V 2<N> 


( 25 ) 


and is typically less than 1*. 

The average surface slope along the satellite ground track can also be estimated by 
differencing the range measured with two different laser pulses (see Fig. 4b). 
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The slope variance using this method is 

Var(tan S x ) = Y^l) + yar(z 2 ) = c_ 2 Var(T D ) _ 
(Ax) 2 2(Ax) 2 

2F 

<N>(Ax) 2 ^ ar + z2tan2 S (tan 2 0x + Var(A<()))] . 


( 27 ) 


This method is potentially very accurate and the error can be as small as 0.1* depending on Ax. 
However, the horizontal separation of the two measurements must be small enough that the 
computed slope represents the true slope trend of the surface. 
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BACKGROUND NOISE, THERMAL 
SAMPLING EFFECTS 


NOISE, QUANTIZATION NOISE 


AND 


The results of the previous sections were derived by neglecting the effects of background 
noise, quantization noise and sampling effects of the A/D converter (see Fig. 1). The A/D 
converter model is illustrated in Fig. 5. A* is the current amplitude corresponding to the kth 
quantization level and At is the sampling interval or range bin length. T is the total observation 

interval. If we neglect speckle noise, the mean and variance of the detector current are 
approximately 

E[i(t)] = G h(t) * [p(t) + p b ] = Ge p(t) + Gept, (2J 

Var[i(t)] = G 2 F h 2 (t) * [p(t) + p b ] + CTlh 2 = [ p(t ) + Pb ] + CTth 2 (2 < 

where w, is the background noise count rate (photons/sec), F is the excess noise factor of the 

detector and ct^ 2 is the thermal noise variance. The mean and variance of the quantized current 
are 


E[i q (mAt)] ~ Ge p q (mAt) + GefXhq 


(30) 



c 

l 

A * quantizer step size 
i Ak - kth quantization level 

At - sampling interval 
T - range computer observation interval 

S 



u 


A 

D 

Cu 

D 

O 

Ak ! 


u 

O 

Fl \ 

o 

<D 

4-* 


T / \ 

D 

Q 


H-At ^ 



-H j ** 

- T H 


Time 


Figure 5. A/D converter model for the laser altimeter. 


Var[i q (mAt)] = [pq(mAt) + M-bq] + <*th + \ 2 ’ (31) 

where the subscript q is used to denote the quantized current and pulse amplitudes and A 2 /12 is 
the quantization noise variance. Note that (31) was derived by assuming the A/D converter 
employs a uniform quantizer, i.e. all the current steps are equal (A = A k +i - A k for every k). 

Equations (30) and (31) can now be used with the previous results to estimate the variances 
of the signal level, propagation delay and pulse width. If we neglect speckle noise, the variances 

are 

Var(N) = F<N> + ^ [F(i b At + (athAt/Ge) 2 + (AAt/Ge) 2 /12] (32) 

Frr „2 T [FqbAt + (a t h At/Ge) 2 + (AAt/Ge) 2 /12] ^ At 2 ,,,, 

Var(Tp) = + At “5? ^ 12 

2FCn 4 T [ F tibAt + (athAt/Ge) 2 + (AAt/Ge) 2 /12] / T 4 T 2 O p 2 4 \ At 4 

V ^>P 2 )=^T + S ^2 U> 6 Pj 144 

(34) 


where a p 2 is given by (18). 
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If Var(Gp) « <jp 2 then 


Var(a p ) = 


Var(q D 2 ) 

4a p 2 


(35) 


and (34) can be rewritten as 


Varfo i _ F ° p2 4 . T [ FtibAt + (cJ th A/Ge) 2 + (AAt/Ge) 2 /12l ( T 4 T 2 On 2> \ 
*" - 2<N> A, <N>* + 


+ 


At 4 

576 a p 2 ' 


(36) 


Notice that the pulse delay and pulse width variances are related to the length of the observation 
period T and the sampling interval At. Because background noise, thermal noise and 
quantization noise are present in the receiver output even when there is no optical pulse, it is 
important to limit the interval over which the pulse parameters are computed to just the period 
when the pulse is present. Equations (33), (34) and (36) were derived by assuming the 
observation interval was (T p - T/2, T p + T/2). The finite range bin width introduces additional 
errors that are related to At. Note that cWGe and A/Ge are the equivalent photon count rates 
associated with thermal noise and quantization noise. 

In a well-designed system the signal level is large compared to background noise so that the 

dominant error source is photon noise in the detected laser pulse. The background count rate for 
solar illumination is given by 

W> = ^ hv ^ S ° (X) ^ m2 tan2 ( 0FOV )] ~2 (37) 

where 

Soft) = solar spectral irradiance measured at the top of the earth ’ s 
atmosphere (Wm~ 2 nm* 1 ) 

AX = receiver optical filter bandwidth (nm) 
t'FOV = telescope half width field-of-view angle (rad) 
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By comparing (37) and (9) and using (32) - (34) it is not difficult to show that background noise 
is negligible provided 

TS 0 (X) AX itr 2 tan 2 (0FOv) « E 

where E is the laser pulse energy. Typically the observation interval extends over the full width 
of the received pulse measured at the e~ 2 points (i.e. T p ± 2o p ) and the receiver field-of-view 
matches the laser divergence angle at the e -2 intensity points so that 


T = 4o p 


(39) 


0FOV = 20t 


(40) 


and (38) becomes 

16 itr 2 0T 2 O p S 0 (X) AX « E . (41) 

In the visible and near infrared region of the spectrum near dawn and dusk S 0 ~ 10 2 
Wm- 2 nm-l. For the case where r = 400 km, 0 T = lmrad, AX = 1 nm and <J p = 10 ns, the left- 
hand side of (41) is less than 1 mj. In this case, background noise will be negligible provided the 
laser pulse energy is much larger than 1 mj. If we use (40) in (37), the background count rate can 

also be expressed as 

p b = T a 2 3 - S 0 (X) AX 0T 2 Ap(<|>) . (42) 

hv 

Now consider a baseline A/D converter design in which the current step size is selected by 
assuming the worse case where all the photons are detected in 1 range bin. 


<N> _ AAt 
2 n Ge 


( 43 ) 
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where n is the number of digitizer bits. We assume the observation interval extends over the full 
width of the received pulse measured at the e~ 2 points (i.e. T p ± 2a p ) so that T = 4 a p . If we 
neglect background noise and thermal noise then the variance of N, T p and a p 2 become 


Var(N) = F<N> + \ ^ 

3 At 2 2n 


(44) 


Var(T p ) = 


Fo p 2 

<N> 


+ 4 ^2 a P 2 . At 2 

9 dt 2 2n 12 


(45) 


Var(c p 2 ) = 


2Fc p 4 

<N> 


. 23 Op Op 4 At 4 
+ 45 ^f^ + T44 • 


(46) 


The A/D converter should be designed so that the errors associated with sampling and 
quantization noise are no larger than the errors associated with photon noise. From Eq. (45) 
we obtain 




~ 2 12 . 


(47) 


(48) 


These equations suggest that the A/D convener should be designed so that the sampling interval 

At is approximately 1/10 the rms width of the received pulse and the number of digitizer bits is 6 
or larger. 
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8. EXPECTED PERFORMANCE CAPABILITIES OF SEVERAL SATELLITE 

LASER ALTIMETER SYSTEMS 

Several laser altimeters are currently under development by NASA for deployment during the 
next decade. The most sophisticated of these is the Geodynamic Laser Ranging System (GLRS). 
GLRS is being developed for the Earth Observing System with a planned launch early in the next 
century. It includes a Q-switched mode-locked Nd:YAG laser operating at the fundamental 
(1064 nm), doubled (532 nm) and tripled (355 nm) frequencies. The altimeter channel utilizes the 
fundamental 1064 nm wavelength and operates with 150 mj/pulse @ 40 pps. The detector is an 
avalanche photodiode. The GLRS altimeter will be used to study ground and ice sheet 
topography, surface roughness and to measure sea state. The major system specifications for the 

GLRS altimeter are listed in Table 4. 

The Lidar In-Space Technology Experiment (LITE) is a shuttle pallet experiment scheduled 
for launch in late 1993 [Couch et al., 1991]. The system includes a flashlamp pumped Nd:YAG 
laser operating at the fundamental, doubled and tripled frequencies. The telescope is a 0.95 m 
diameter beryllium mirror. The 1064 nm altimeter channel was designed to study the reflectivity 
of clouds, land and the ocean. The laser pulse energy is 486 mj @ 10 pps. The range resolution 
is a rather modest 15 m which is limited by the 10 MHz waveform digitizer. The major system 
specifications for LITE are also listed in Table 4. 
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Table 4 NASA Spaceborne Laser Altimeters 


Instrument 

LITE 

TMLA 

GLRS 

MOLA 

LOLA 

Planet 

Earth 

Earth 

Earth 

Mars 

Moon 

NASA Mission 

Shuttle 

Earth Probe 

EOS 

Mars Observer 

Lunar Observer 

Project Start Date 

1988 

1994 

1989 

1989 

1994 

Launch Date 

1993 

1998 

2000 

1992 

1998 

Mission Duration 

7 days 

3 yrs 

5 yrs 

2 yrs 

2 yrs 
175 km 

Orbital Altitude 

300 km 

400 km 

705 km 

400 km 

Laser Footprint Dia. 

300 m 

200 m 

100 m 

160 m 

30 m 

Laser Parameters 




Wavelength 

1064 nm 

1047 nm 

1064 nm 

1064 nm 

1064 nm 

Energy/pulse 

486 mj 

15 mj 

80 mj 

40 mj 

10 mj 
40 pps 

Pulse rate 
Pulse width 

10 pps 

1050 pps 

40 pps 

10 pps 

(FWHM) 

27 ns 

5 ns 

5 ns 

8 ns 

5 ns 

Divergence 0? 




HW@e“^ 

250prad 

125 prad 

35.5 (trad 

100 prad 

43 prad 

Telescope Diameter 
Obscuration 

0.95 m 

0.9 m 

0.6 m 

0.5 m 

0.5 m 

Diameter 

0.31 m 

0.25 m 

0.15 m 

0.12 m 

0.12 m 

Pointing Control 

20 mrad 

50prad 

50 prad 

3 mrad 

3 mrad 

Pointing Jitter 

3 mrad 

10 prad 

10 prad 

1 mrad 

20 prad 

Range Resolution + 
Overall Optical 

15 m 

0.1 m 

0.1 m 

1.5 m 

0.2 m 

Efficiency* 

10% 

20% 

20% 

20% 

20% 

Received Pulse 
Statistics 




< N > 

87,000 

2,800 

2,200 

4,700 

6,200 

K, 

4.4 x 10 s 

1.1 x 10 5 

3,700 

2.1 x 10 4 

3,800 

SNR 

153 

28.2 

23.2 

35.5 

34.8 


+ Range bin length Az = cAt/2 
♦Includes APD quantum efficiency (-50%) 


The Mars Observer Laser Altimeter (MOLA) is currently in the final testing phase and will 
be launched in 1992. MOLA is designed primarily to map the topography of Mars but will also 
provide some information on surface roughness. The altimeter includes a 40 mj/pulse, 10 pps 
NdrYAG operating at 1064 nm and a 0.5 m diameter telescope. LOLA, the Lunar Observer 
Laser Altimeter, is essentially an upgraded version of MOLA which is planned for launch in 

1998 to map the lunar topography. The system specifications for both MOLA and LOLA are 
listed in Table 4. 
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The Topographic Mapping Laser Altimeter (TMLA) is designed specifically to measure land 
and ice surface elevations globally. TMLA is envisioned as a small dedicated satellite of the 
Earth Probe class and is planned for launch in a polar sun-synchronous orbit in 1998. To obtain 
high resolution coverage, the Nd:YLF laser array will operate at 15 mj/pulse @ 1050 pps. The 
pulse length is 5 ns. The TMLA system specifications are summarized in Table 4. 

To evaluate the performances of these altimeters it is necessary to make some assumptions 
about the target characteristics. For simplicity we assume a nominal value of 0.3 for the diffuse 
reflectivity of the earth, moon and Mars. At zenith for excellent visibility, the 1-way 
transmission through the earth’s atmosphere is approximately 0.95 at 1064 nm. To be 
conservative we assume the 1-way atmospheric transmission is 0.7 for the earth and 1.0 for the 
moon and Mars. The mean detected photon count per pulse is listed in Table 4 for each 
altimeter. Also listed is the signal-to-noise ratio (SNR) which is defined as 

SNR = <N>/Std(N) = (49) 

The excess noise factor of the APD was assumed to be 3.5 for each system. LITE has the highest 
SNR, primarily because the large pulse energy and telescope provide a very high signal level and 

the large laser divergence minimizes speckle noise. 

Ranging accuracy is influenced considerably by the terrain slope and surface roughness. 
Representative terrain statistics for the earth’s surface are listed in Table 5. These data were 
obtained by Harding et al. [1991] by analyzing digital terrain elevation data with 90 m spatial 
resolution. The slope values correspond to the mean north-south and east-west surface slopes. 
The £ 50% and L 90% entries refer to the surface slope and RMS roughness values representing 
the 50% and 90% cumulative distribution levels. For example, 50% of the north-south and east- 
west surface slopes in high relief areas were less than 13.6' while 90% were less than 28.2 . 
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Table 5 Representative Terrain Statistics 


Relief 

Geom orphic Features 

Surface Slope 

RMS Roughness 



Z50% 

Z90% 

150% 

190% 

Low 

Glaciated Continental Shield 

mmm 

f 1.6' 

0.5 m 

0 8m 


Shallowly Incised Drainage 


mgm 

0.9 m 

1.9 m 


Ice Cap 

0.90* 


0.3 m 

1 ,0 m 


Moderately Incised Drainage 

0.95* 

Km 

1.4 m 

3.1 m 


Hnean 

0.81° 

mm 

■HEEl 

1.7 m 

Moderate 

Quaternary Volcanic Complexes 

1.0* 



5 8m 


Normal-Faulted Rift Graben 

i.r 

mm 

1.2 m 

4.7 m 


Deeply Incised Drainage 

4.8* 

MM 

5.3 m 

14.1 m 


Fold and Thrust Orogenic Belt 

3.2* 

17.0* 

4.0 m 

11.1m 


i mean 

2.5' 

12.5' 

3.0 m 


Hl 'g h 

Convergent Mountain Front 

En 

msm 

6.3 m 



Table 6 Altimeter Performance at Nadir for Low Relief Terrain 
Surface Slope = 0.8° RMS Roughess = 0.8 m 


Instrument 

LITE 

TMLA 

GLRS 

MOLA 

LOLA 

RMS Received Pulse Length 






System Effects (cm) 

470 

32 

32 

67 

32 

Beam Curvature (cm) 

2 

! 0.6 

0.09 

0.4 

0 03 

Surface Roughness (cm) 

80 

80 

80 

80 

80 

Surface Slope (cm) 

150 

99 

49 

79 

15 

Total (cm) 

500 

130 

100 

130 

90 

RMS Range Error 






System Effects (cm) 

430 

3.1 

3.2 

43 


Pointing Jitter (cm) 

1780 

7.9 

14 

790 

Surface Roughness (cm) 

0.52 

2.8 

3.5 

2.3 


Surface Slope (cm) 

0.95 

3.5 

2.0 

2.2 

Total (cm) 

1830 

9.6 

15 

790 

9.3 

RMS Pulse Length Error 






System Effects (cm) 

188 

0.8 

0.9 

73 

0.57 
I 3 

Surface Roughness (cm) 

0.36 

2.0 

23 

1.5 

Surface Slope (cm) 

0.66 

2.5 

1.4 

L5 

0.25 

Total (cm) 

188 

3.3 

2.8 

7.6 

1.4 
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Table 7 Altimeter Performance at Nadir for Medium Relief Terrain 
Surface Slope = 2.5° RMS Roughess = 3.0 m 


Instrument 


RMS Received Pulse Length 

System Effects (cm) 

Beam Curvature (cm) 
Surface Roughness (cm) 
Surface Slope (cm) 

Total (cm) 


RMS Range Errors 

System Effects (cm) 
Pointing Jitter (cm) 
Surface Roughness (cm) 
Surface Slope (cm) 

Total (cm) 


RMS Pulse Length Errors 

System Effects (cm) 
Surface Roughness (cm) 
Surface Slope (cm) 

Total (cm) 



Table 8 Altimeter Performance at Nadir for High Relief Terrain 
Surface Slope = 13.6° RMS Roughess = 6.3 m 


Instrument 


RMS Received Pulse Length 

System Effects (cm) 

Beam Curvature (cm) 
Surface Roughness (cm) 
Surface Slope (cm) 

Total (cm) 


RMS Range Errors 
System Effects (cm) 
Pointing Jitter (cm) 
Surface Roughness (cm) 
Surface Slope (cm) 

Total (cm) 


RMS Pulse Length Errors 

System Effects (cm) 
Surface Roughness (cm) 
Surface Slope (cm) 

Total (cm) 
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The X 50% statistics in Table 5 were used to compute the RMS received pulse lengths 
c/2 [E(o p 2 )] 1/2 and the single-shot range errors c/2 Std(T p ) and pulse length errors c/2 Std(a p ) 
for the LITE, TMLA, GLRS, MOLA and LOLA altimeters. The results for low, medium and 
high relief terrain are tabulated respectively in Tables 6, 7 and 8. The simplified formulas used 
to compute the data in Tables 6-8 are listed in Table 9. 

The data in Tables 6-8 show that pulse broadening is dominated primarily by terrain effects. 
For low and medium relief terrain, surface roughness effects are the most significant for GLRS, 
MOLA and LOLA while surface slope effects are most significant for TMLA because of its 
larger divergence angle. System effects, viz. the relatively slow 10 MHz waveform digitizer, 
make the largest contribution to pulse broadening for the LITE altimeter over low and medium 
relief terrain. For high relief terrain, slope effects make the largest contribution to pulse 
broadening for the LITE, TMLA, GLRS and MOLA aldmeters. For LOLA, slope effects are not 
as important as roughness effects because the orbital altitude is low and the laser divergence is 
small. 

Because the signal levels for all five altimeters are quite high, pulse broadening contributions 
to the range error are relatively small and typically only a few cm. The major source of range 
error is pointing jitter. Pointing jitter range errors are most severe over high relief terrain where 
the slopes are large. Because pointing jitter is large for LITE (3 mrad) and MOLA (1 mrad), the 
single-shot range error for these systems varies between about 1 and 30 m depending on the 
terrain statistics. The single-shot range error for TMLA, GLRS and LOLA are only 10-50 cm 
over low and medium relief terrain because the pointing stability for these altimeters is 100 |i.rad 
or less. Terrain effects are the major source of errors in the single-shot pulse width 
measurements for TMLA, GLRS, MOLA and LOLA while the waveform digitizer resolution is 
the major source of pulse width errors for LITE. 
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Table 9 Simplified Formulas Used to Compute 
the Entries in Tables 6, 7 and 8 


RMS Received Pulse Length c/2 [E(a p 2 )] 1/2 

System Effects j 

c/2 (at 2 + At 2 /12) 1/2 

Beam Curvature 

z tan 2 9r 

Surface Roughness 

Std(A^) 

Surface Slope 

V 2 z tanBj tan S 

RMS Range Error c/2 Std(T p ) 

System Effects 

c/2 (at 2 + At 2 /12) + At 2 /12] 1/2 

Pointing Jitter 

V 2 z tan S Std(A<|>) 

Surface Roughness 


Surface Slope 

(^*2K;) l ' 2,,iz “ er ' anS 

RMS Pulse Length Error c/2 Std(<T p ) 

System Effects 

c/2 ^2<N> + At2/l2 ) + 576 E(a p 2 )l 

Surface Roughness 


Surface Slope 

( dk ) 1/2 z tan0T tanS 


9. CONCLUSIONS 

Satellite laser altimeters are capable of making important contributions to our knowledge of 
the topography of the earth and planets. Surface elevations, slopes and roughness can be 
measured with high precision and accuracy. Well designed systems are capable of achieving 
single-shot ranging accuracies of a few centimeters to a few meters depending on the terrain 
characteristics. Excellent signal levels can be obtained using diode-pumped solid state lasers 
such as Nd:YAG with pulse energies of a few tens of mj and telescope diameters of a few tens of 
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cm. Ranging accuracy is critically dependent on the pointing accuracy and stability of the 
altimeter especially over high relief terrain where the surface slopes are large. At typical orbital 
altitudes of several hundred km, single-shot range accuracies of a few cm can be achieved only 
when the pointing jitter is on the order of 10 (irad or less. Of course ranging accuracy can be 
improved through pulse averaging. 

For earth observations, perhaps the most important application for laser altimeters is the 
measurement of the Greenland and Antarctic ice sheets. Knowledge of the seasonal and long 
term changes in ice sheet volume can play a significant role in testing and validating global 
climate models. Ranging accuracies of a few tens of cm or better are required for this 
application. Over the oceans, laser altimeters can also be used to measure sea state and surface 
wind speeds, while accurate measurements of surface roughness and topography over land have 
numerous applications in the geological sciences. 
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